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Abstract 

We investigate theoretically the transport of non-interacting electrons 
through an Aharanov-Bohm (AB) interferometer with two quantum dots 
(QD) embedded into its arms. In the Coulomb-blockade regime, trans- 
port through each QD proceeds via a single resonance. The resonances 
are coupled through the arms of the AB device but may also be coupled 
directly. In the framework of the Landauer-Biittiker approach, we present 
expressions for the scattering matrix which depend explicitly on the en- 
ergies of the two resonances and on the AB phase. We pay particular 
attention to the crossing of the two resonances. 

PACS numbers: 72.10.Bg, 72.20.Dp, 72. 20. My, 72.15.Qm 

1 Introduction 

The crossing of two Coulomb-blockade resonances was studied in two recent 
experiments ^0- In both cases, two quantum dots (QD) were imbedded into 
the arms of an Aharanov-Bohm (AB) interferometer. By changing the param- 
eters of the experiment (various gate voltages and the magnetic flux through 
the AB device), it was possible to study the crossing properties of two iso- 
lated Coulomb-blockade resonances, one each due to one of the two QDs. In 
the present paper we present a theoretical framework for the analysis of both 
experiments. 

Figure 1 shows a schematic representation of both experiments. The AB 
ring contains the two QDs labelled QDL and QDR where L and R stand for 
left and right, respectively. The QDs are separated by barriers from the rest of 
the AB device. The latter consists of two parts. In Figure 1, the lower (upper) 
part is labelled 1 (2, respectively). Both parts are coupled to the outside world 
by a number of leads. In Figure 1, this number is two (three) for part 1 (part 
2, respectively). In our theoretical treatment, the number of leads coupled to 
each part will be arbitrary. Typically, one of the leads coupled to part 1 (part 2) 
serves as source (sink, respectively) for the electrons. While the two QDs are not 
coupled directly to each other in the first experiment p^, such a coupling does 
exist in the second experiment This coupling is indicated schematically by 
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the dotted horizontal line representing the wire connecting QDL and QDR. In 
Ref. 0, the strength of that coupling was controlled by a further gate. Figure 1 
does not show the plunger gates which make it possible to control the energies 
of the Coulomb blockade resonances in either QD. Thereby it is possible to have 
the energies of both Coulomb blockade resonances coincide. Experimentally, 
such crossings are seen in three-dimensional plots of the conductance versus 
the plunger gate voltages Vl and Vr applied on QDL and QDR, respectively. 
Each Coulomb-blockade resonance corresponds to a ridge. The ridges of res- 
onances in QDL (QDR) run essentially parallel to Vr (Vl, respectively). The 
crossing of two such ridges marks the crossing of two Coulomb-blockade reso- 
nances. The coincidence of two resonances also affects the interference pattern 
of the transmission of an electron through the AB device. This pattern depends 
upon the magnetic flux $ through the device. The flux is due to a homogeneous 
magnetic field perpendicular to the plane of the drawing. We are interested in 
weak magnetic fields only. (We recall that for a complete AB oscillation, the 
magnetic field strength typically changes at most by several ten mT.) There- 
fore, we take into account only the AB phase due to the magnetic flux and 
neglect the influence of the magnetic field on the orbital motion of the electron. 
Gauge invariance then allows us to link the AB phase to the passage of the 
electron through a particular part of the AB device. In the absence of a direct 
coupling between the two QDs (i.e., without the dotted line in Figure 1), we 
choose the barrier separating QDL from part 1. Whenever the electron leaves 
(enters) QDL for part 1 (from part 1), it picks up the phase factor exp(2«7r<I>/<I>o) 
(exp(— 2z7r<I>/<I>o), respectively) where $0 is the elementary flux quantum. For 
brevity, we write the phase factor as exp(i(f>). In the presence of a direct link 
between the two QDs, the topology of the AB interferometer changes from that 
of a ring to that of a figure eight, and we use a different convention in Section|3| 
In Section [5] we define the Hamiltonian for the system. In Section [3] we use 
the Landauer-Biittiker approach and present the generic form of the scatter- 
ing matrix which describes the experimental setup of Ref. . In Section 21 
this scattering matrix is analysed especially with regard to the crossing of two 
Coulomb-blockade resonances. In Section [5] we generalize our treatment to in- 
clude the setup of Ref. [2]. In Sectional we list the approximations and summa- 
rize our approach and results. Moreover, we address some of the approximations 
made. In particular, we discuss the neglect of the mutual Coulomb interaction 
between the two electrons which are added to the system as the resonances 
become populated, and that of the Coulomb interaction between each of these 
electrons and those on the dots. We also address the role of the spins of both 
quantum dots and of the two added electrons. Throughout the paper, we disre- 
gard temperature averaging for simplicity. Likewise, we disregard decoherence 
effects although these are known to play some role in the actual experiments. 
We do so because part of the transport through the device is known to pro- 
ceed coherently. Only this part will display a dependence on the AB phase. 
Moreover, decoherence has been thoroughly discussed in the literature, see, for 
instance, Ref. 
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2 Hamiltonian 

In defining the Hamiltonian of the system, we proceed in full analogy to Refs. |3J 
[5]. These papers addressed the AB phase for a single QD placed in one of the 
arms of an AB interferometer. (For a review of work on this problem, sec 
Ref. [S])- We introduce fictitious barriers separating parts 1 and 2 of the AB 
device from the attached leads. Likewise, we consider parts 1 and 2 as separated 
from the two QDs. We impose boundary conditions on all these barriers such 
that as a result, we obtain self-adjoint single-particle Hamiltonians i?i ca d for the 
leads, H\ and H2 for the now separated parts 1 and 2, and H^ and Hr for the two 
QDs labelled QDL and QDR, respectively. Here H\ ea< i possesses a continuous 
spectrum while the spectra of Hi, H2, £2l and £Zr are discrete. We label the 
leads attached to part 1 (part 2) by s = 1, . . . , S (by t = 1, . . . , T, respectively). 
The transverse modes (channels) in lead s (t) are labelled a — 1, . . . , N s (a — 
1, . . . ,Nt, respectively), and correspondingly for the creation and annihilation 
operators and c. The associated energies are labelled e. The eigenvalues of 
Hi (H2) are labelled Eij (E2j), with j = 1, . . . , 00 and associated creation and 
annihilation operators Cy (cL)> and cij (c2j, respectively). We assume that 
transport through cither QD occurs in the Coulomb-blockade regime where the 
intrinsic widths of individual resonances are small compared to their spacings. 
(The spacing includes, of course, the charging energy). We also assume that the 
temperature is small in comparison with the spacings. Under these conditions, 
it is legitimate to assume that transport through either QD is dominated by a 
single Coulomb-blockade resonance. We believe that this situation is met or 
nearly met in the experiments of Refs. 012- Thus, we admit only a single 
bound state with energy El (El) in QDL (QDR), with associated creation and 
annihilation operators d L (d^) and dh (g?r, respectively). The energies E\, and 
El include the charging energies. Altogether, we have 
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flu, = ^{Vij^cljd-L+h.c.) , 
j 

i 

ffm = ^(V^-.flcJjdR + h.c.) , 

#2R = X^« C 2A + h - C ") • ( 2 ) 

j 

The direct coupling of QDL and QDR (dashed line in Figure 1) is given by 

H hR = V hK (4dK + h.c.) . (3) 

In the absence of any direct coupling between QDL and QDR (Vlr = 0) we use 
gauge invariance to put the entire AB phase onto a single one of the barriers. 
Without loss of generality we choose the barrier separating QDL and part 1. 
Then, all the matrices V in Eqs. J5J are real and symmetric except for Vij-i, 
which obeys 

Vij-L eKp(-i(f>) = V L -ij exp(i(/>) = vi j]L (4) 

with vij-L real and symmetric. For Vlr 7^ 0, a modification is neccessary and 
discussed in Section below. The Hamiltonian H of the system is the sum of 
the terms defined by Eqs. Q to ©. We have not considered the possibility of 
spin-orbit coupling on either QD. 

We have been very explicit in the construction of H. The reason is that 
we wanted to show that H is a sum of single-particle Hamiltonians. This fact 
allows us to use the Landauer-Biittiker approach to describe transport through 
the system. The ensuing use of the scattering matrix enables us to display 
explicitly the phase- and energy-dependence of the conductance coefficients. 
We have omitted the spins of as well as any possible interaction between the 
two electrons which will eventually populate the two resonances caused by El, 
and Sr. These points are taken up in Sectional 

3 Scattering Matrix: Ring Topology 

The transport through the device is described by the Landauer-Biittiker formula 
I s = ^G ss ,Vv +Y J GstV t , s= 1,...,S , 

s' t 

h - J2 Gtt,Vt ' +J2 GtsV * > * = 1 =---> T ' (5) 

t' s 

Here I s (It) is the current through lead s (lead t), respectively, and V s (Vt) is the 
voltage applied to that lead. The conductance coefficients G ss > are given by 

N s N a , 

G <* = 2 E l\Ssa;s>a'(E : $)| 2 - 5 SS ,} (6) 
a— 1 a' — l 
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and correspondingly for the index combinations (st), (ts) and (tf ). The sym- 
bol S S a;s'a'(E, <&) denotes the element of the scattering matrix S(E,§) which 
connects channel a in lead s with channel a' in lead s' at energy E and mag- 
netic flux <!>. Time-reversal symmetry requires the scattering matrix to obey 
the relation 

S(E,$) = S T (E,-<f>) (7) 

where T denotes the transpose. 

The observable O describing a given experimental setup is determined by 
the experimental arrangement chosen (which of the leads are grounded, and in 
which of the leads is a current measured). For any such setup, O will be given as 
a rational function of the conductance coefficients G. This follows directly from 
Eq. 0. In order to present a general framework useful for the analysis of any 
such experiment, we focus attention on the scattering matrix S. With the help 
of the formulas for S given below, it is possible to work out the dependence of the 
G's and, hence, of O on the energies of the two Coulomb-blockade resonances, 
and on the AB phase <&. 

It is possible to derive the form of S from the Hamiltonian H . This can be 
done along the lines of Refs. 0J[7|. We do not follow this course here because the 
explicit solution involves some lengthy algebra. Rather, we simply present the 
result which we believe to be intuitively obvious. In this and the next Section, 
we focus attention on the ring geometry and put Vlr = 0. 

The scattering matris S can be written as the product of three unitary 
matrices, 

S(E, $) = US^ (E, <1>)U T . (8) 

Without any coupling between each of the QDs and parts 1 and 2 of the AB 
device (this condition can be met experimentally by increasing the heights of 
the two barriers defining each QD), the resonant part S^ IQS \E, is equal to 
the unit matrix, and S(E, $) is, thus, equal to UU T . The form of the latter 
matrix follows from the observation that parts 1 and 2 are unlinked. A uni- 
tary scattering matrix (S^) describes the non-resonant electron transport 
through unlinked part 1 (unlinked part 2, respectively). We assume that the 
energy dependence of both matrices is smooth over the energy interval defined 
by the widths of the two Coulomb-blockade resonances introduced below. We 
accordingly neglect the energy dependence of both and ifv 2 ' . Moreover, 
both matrices do not dependent on the magnetic flux <!>, see the remark at the 
end of Section H] Time-reversal invariance then implies that both and S^ 2%> 
are symmetric. Thus, we can write for i = 1,2 

=U®[U®] T . (9) 

Eq. JJjJ holds for every unitary and symmetric matrix. As explained in Refs. [HI 
EJ, the unitary transformation {/W accomplishes the transformation from the 
space of physical channels to the space of eigenchannels. We accordingly write 
the matrices explicitly in the form UsV-a and uf£!g. Here UsV-.a is the 
product of an orthogonal matrix O^J- a which diagonalizes the symmetric matrix 
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and of a diagonal matrix with entries exp(i8a ) where the 8& s are the 
eigenphaseshifts of and similarly for S^ 2 \ The index a (/3) runs from 1 to 
N\ (to A^2, respectively). Here the total number of channels N± in part 1 (N2 
in part 2) is given by N\ = J2 S (by N 2 = J2t Nt, respectively). The matrix 
U is defined in the total space of N = N% + N2 channels. It is block-diagonal 
and given by 

u ={ U o ik) ■ ( 10 ) 

Inspection of Eq. (JSJ) shows that for S^ Tes > = In, the unit matrix in N dimen- 
sions, the scattering matrix S is block-diagonal and consists of the two matrices 
S* (1) and S" (2) , as it must. 

It is now obvious that S'^ rcs - ) differs from the unit matrix by terms which 
represent the two Coulomb-blockade resonances, one each in QDL and QDR. 
Moreover, it is also clear that 5( ros ) is defined in the space of eigenchanncls of 
both and S<- 2 \ In this space, the coupling matrix elements W p p describing 
the hopping of an electron from the resonance in QDP (with P = L or R) to the 
eigenchannel p (with p = a for part 1 and p = j3 for part 2) can be shown [H] to 
be real, save for the AB phase. We accordingly have for p = 1, . . . , N and P = 
L,R 

Wp p = Wp p = W p p unless P = L and p = a (11) 

while 

Wh a exp(z0) = WaL exp(-z^) = w Lce (12) 

with wl« real. We note that the Wp p 's differ from but are linear in the Vp ;so 's 
and Vp ; t a 's introduced in Section |2| 

We can now express S^ ^ in terms of the matrix elements Wp p , and of the 
energies Ep of the two Coulomb-blockade resonances. The latter can be varied 
experimentally by changing the plunger gate voltage on either QD. We observe 
that the matrices Wp p map the space of N eigenchannels onto the space of the 
two Coulomb-blockade resonances, and vice versa for W p p. The matrix 5^ res -' 
takes the form 

Sffi = 6 pp' ~ 2 ™J2 W pF [D- l } FP ,W P , p , . (13) 
pp' 

The two-by-two matrix -Dpp' has the form (P = L, R) 

D P p, = 5p P , [E - E P ] + in Wp p W p p, . (14) 

p 

Eqs. H6I8I13I14JI constitute the central result of this Section. It is easy to check 
that the scattering matrix defined by these equations is unitary and obeys 
Eq. Q. 

For the benefit of the reader, we rewrite the S'-matrix in a form which 
displays more clearly the physical role of the matrices J/W with i = 1,2. We 
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define the complex coupling matrix elements 



W ta ., P = ^^>p. (15) 



Then S takes the form 
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The matrix D has the same form as in Eq. I|14fl but can also be written as 

£>P P / = Spp, [E - Ep] + in Wp W S a-P' + ITT J2Wp;taWta;P> ■ (17) 

sa ta 

The transformation (|15|l introduces complex matrix elements W which guaran- 
tee unitarity of S in the presence of the non-diagonal unitary matrices and 

4 Analysis. Crossing of Two Resonances 

The effect of the two resonances which dominate the scattering matrix is con- 
tained entirely in the matrix D defined in Eq. (|14|l . It is useful to display D in 
matrix form, 

n _ (e - e l + (</2)r L (z/2)r LR \ 

I i; '-nr. . r . r . i , ; .'->'ir . J • 

where 
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Our explicit notation combined with Eqs. and (|12fl shows that IY and Tr 
are real, positive and independent of the magnetic flux $, and that the only 
dependence on $ occurs in IYr and in Trl- The latter two quantities are 
complex and related by 

Plr — Trl • (20) 

We use Eqs. ijTHj) to lj2"U)l to display the structure of certain elements of the 
scattering matrix S. We recall that S decays into two independent scattering 
matrices and whenever we have Wp p = for all P, p. Parts 1 and 2 of 
the AB interferometer are linked only by the two Coulomb-blockade resonances 
with energies E^ and E R . If the two resonance energies are sufficiently different 
so that 

\E h -E n \ »r L ,r R) |r LR | , (21) 

we can use perturbation theory in Tlr to invert D. Keeping only the lowest- 
order terms in the expansion, we find 

s tb;sa = u & [ W ^\ E ~ + (*/2)r L ]- 1 ^ LQ 

pa ^ 

+Wp R [E -E R + (i/2)r R ]" 1 W Ra ^ Uj$ a . (22) 

The two amplitudes on the right-hand side of Eq. i|22|) can be interpreted in 
terms of two paths of the electron on its way from part 1 to part 2. The electron 
may pass either through QDL (first term) or QDR (second term). As it passes 
through QDL, it picks up the AB phase contained in Wt, a - This phase will affect 
the interference pattern due to the product of the amplitudes corresponding to 
the two paths. Whenever inequality l|21|) holds, the electron will not complete 
one or several loops within the AB ring as it passes from part 1 to part 2. It is 
instructive to consider also the terms of next order. These terms are given by 

(w pl [e -e l + ^^rLr^LRji; - e r + ^/2)r R ]- 1 ^R Q 

0a ^ 

+W (3K [E -E R + ( l /2)r R ]- 1 rR L [s -E L + {i^T^WLc^U^l . 

(23) 

The path associated with the first amplitude leads the electron first through 
QDR and then through QDL, and vice versa for the second amplitude. With 
Plr given by Eq. (|19|l . we see that along the first path QDL can be reached 
from QDR cither via part 1 or via part 2, and correspondingly for path 2. In the 
first (second) case, the AB phase does not (does) contribute to the scattering 
amplitude. This is correct because only in the second case does the electron 
complete a loop around the AB ring. A similar analysis of S sa - S 'a' and of 
Stb-t'b' shows that whenever inequality (|21|l holds, the scattering is dominated 
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by and S^ 2 \ respectively. The consecutive passage through both Coulomb- 
blockade resonances is strongly inhibited. We conclude that inequality 
defines a fairly uninteresting regime of parameters of the problem. 

Interest, therefore, focusses on the regime where this inequality does not hold 
and where our perturbation expansion is not appropriate. This is the regime 
where the two Coulomb blockade resonances may cross. We shall see that the 
crossing displays novel features. Prior to calculating the exact result, it is useful 
to visualize the outcome in terms of a perturbation expansion in powers of Tlr 
and Trl ■ This expansion generates terms of the same form as in formula (I23[) but 
of higher order in Tlr and Trl- Each propagator \E— Ep + (i/2)Tp]~ 1 occurring 
in the expansion signals a visit of the associated path to QDP with P = L, R. 
The intermittent factors Tlr and Trl signal passage of the electron from QDR 
to QDL and vice versa. The passage may proceed via part 1 or part 2. Thus, 
the perturbation series stands for the infinite number of possibilities to connect 
the channels defined by the indices of the scattering matrix, by paths. These 
paths may loop around the AB ring a number of times, then change direction, 
loop again, change direction back etc. until the electron leaves the AB ring. 
The AB phase picked up by the electron is the sum of all such phases picked 
up in the individual loops and given in terms of the total number of completed 
counter-clockwise loops minus the total number of completed clockwise loops. 
We are about to calculate the form of the scattering matrix by diagonalizing 
the matrix D. This procedure amounts to summing over all the paths just 
mentioned. This is why the AB phase will show up in the denominator of the 
result, see Eq. (|27(l . (Experimentally, decoherence will actually limit the number 
of loops that contribute significantly to the amplitude, see Ref. £Q). 

We simplify the algebra by considering an AB ring which contains two per- 
fectly identical QDs and which itself is perfectly symmetric about a vertical axis 
through the middle of Figure 1. Then, Tl = Tr = T (this defines the width 
P) . We write the complex eigenvalues of the matrix D in the form E — Si with 
i = 1,2. Then 



£i, 2 = ^(E L + E R ~ %T) ± - V(El - £r) 2 - |r LR | 2 ■ (24) 

Let us suppose that we change the resonance energies of both dots in such a 
way that (I?l + Ep) is kept fixed while u = |£r — £r| decreases monotonically 
from an initially large value (in the sense of the inequality $'21\i ). Then, the 
difference \ei — £2! also decreases monotonically. Both resonances approach 
each other, retaining equal widths. The difference of resonance energies vanishes 
when (_El — Ep) 2 = |Tlr| 2 : The two resonances coincide in energy and width. 
We deal with an exceptional point in the sense of Ref. [S]- At this point, the 
system possesses only a single eigenfunction. As we decrease u further, the 
two resonances separate, retain equal resonance energies but acquire different 
widths. At u = 0, the widths differ by ITlrI, the maximum amount possible. 

The value of |Flr| determines both, the value of u where the resonances 
coincide and the maximum difference of their widths. This value depends upon 
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the AB phase $. Indeed, from Eq. (|19fl we have 




47r 2 [(J] w La w aR ) 2 + W L0 W 0R ) 2 





(25) 



The value of |Tlr| 2 oscillates periodically with magnetic flux <I> between the 
maximum value 47r 2 (^ Q WLaWaR + Ye Wl/?W/9r) and the minimum value 
47r 2 (^ Q WhaWaR — Yp Wl^W^r) 2 . From Schwarz's inequality we conclude 
that the widths of the two resonances are always positive. To estimate the 
relative size of T and of |Tlr|, we note from Eq. I|19|) that r = IY = Tr is a 
sum of squares while IYr is a sum over terms which, aside from the AB phase, 
may have either sign. We expect that due to impurity scattering in parts 1 and 
2 of the AB device, the Wp^'s are Gaussian random variables, see Ref. 5 . As a 
consequence, we have T cx TV while ITlrJ fluctuates strongly with a root-mean- 
square variance which grows like \/~N. Thus, the maximum difference of the 
widths of the two resonances is expected to be of the order of T/y/W. 

The possibility of complete coalescence of two resonances displayed above is 
a phenomenon which is opposite to the well-known Wigner-von Neumann level 
repulsion effect for bound states. The latter occurs whenever two bound states 
interact via a Hermitean interaction. Eq. (|18|l shows that in the present case, we 
deal with resonances with complex resonance energies to begin with, and with 
a coupling that is due to a Hermitean interaction multiplied by i, the imaginary 
unit. This unusual form of interaction occurs because the two resonances are 
not coupled directly but via the open channels in parts 1 and 2. Both differ- 
ences contribute towards a behavior which differs from standard Wigner-von 
Neumann level repulsion. Such behavior has been discussed previously in the 
literature. To the best of our knowledge, the coupling of two resonances was 
first studied explicitly by von Brentano et al. J0| m the context of Nuclear 
Physics. This work was followed by an experimental investigation Related 
work was published in Ref. ^2 • Recent work ^2] has focussed on the properties 
of exceptional points. 

To display the features of the exceptional point where (E R — -El) 2 = |Tlr| 2 
and where IY = Tr, we consider two slightly asymmetric QDs for which the two 
resonance widths IY and Tr are not exactly equal. Then, the eigenvalues 
will never coincide exactly. This is seen from the expression of the discriminant 
which now has the value 



Imagine now a change of the the parameters of the system in such a way that 
the argument of the square root describes a closed loop in the complex plane 
around the exceptional point. This could be achieved as follows. We put E R = 
-E h = |r LR | +a/2, T L = -T R = with a, {3 real an d \a\ < |r LR |,|/ 3| < |r LR |. 
The discriminant becomes approximately equal to -\/2|rLR|(a — i(3). Changing 



y/(E h -e r - (i/2)(r L - r R )) 2 - |r LR | 2 . 



(26) 
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a from a small negative to a small positive value while keeping (3 > fixed and 
small, then keeping a fixed and changing [3 from its small positive value to a 
small negative one, then keeping f3 fixed and changing a back to its original value 
and doing, finally, the same for /3 yields a rectangle in the complex plane with 
the exceptional point in its interior. While under this operation the phase of 
(a — i(3) changes by 2-zr, the phase of ^/2|rLR,|(a — i(3) changes only by ir: Under 
this operation, the two eigenvalues £1,2 are interchanged, and so are the two 
eigenfunctions, including an additional phase factor In comparison with 

the work of Refs. ^3, the present system seems to offer an additional degree 
of freedom in terms of the AB phase. The latter determines the value of ITlrI, 
the location of the exceptional point, and the form of the two eigenfunctions as 
linear combinations of the two QD states. 

Unfortunately, all these appealing features have no bearing on the properties 
of the scattering matrix S. This is because the energy E is always real. As a 
consequence, we can never reach the exceptional point, and the determinant 
of the matrix D never vanishes for real values of E. The two eigenfunctions 
of D remain distinct. In view of the recent interest in exceptional points, we 
have nevertheless felt that a discussion of this topic is appropriate in the present 
context. 

We return to the symmetric case. The matrix D can be diagonalized by a 
matrix A so that D — A (EI2 — e)A where e denotes the diagonal matrix 
diag(ei,£2). Using this form in Eq. we obtain 

4p' S) = <W - 2i7T W pP A pk i E - ep 1 ]~ 1 Mp'Wp' P > ■ (27) 
ppiP' 

Inserting this matrix into Eqs. © yields the conductance coefficients and, hence, 
the dependence of any observable on the AB phase. The AB phase appears 
explicitly not only in the eigenvalues £1,2 but also in the matrix A and, of 
course, in some of the W p p's. The matrix A can easily be calculated. Details 
are not given here. 

For the sake of completeness, we discuss the limitations of a two-lead experi- 
ment. These limitations have played a role in previous studies of AB devices jH]. 
We recall that the scattering matrix S is unitary and obeys S T (E, — $) = 
S(E,$). It follows that in general, we have S sa - S ' a ' (E, -$) = S s > a > ;sa (E, $) 
and correspondingly for the lead indices (s, t) and (t,t'). The cases where the 
two lead indices coincide are special and yield S sa;sa < (E, — $) = S sa ' ;sa (E, $) 
and Sta-,ta'{E, — $) = Sta'-.ta{E, $). For the conductance coefficients, this means 
that G^(E,-$) = G w (E,-$) = G vt {E,$) and G st (E, -$) = 

G ts (E, $) while G SS (E, -$) = G SS (E, $) and G tt (E, -$) = G tt (E, $). Unitar- 
ity then shows that for a two-lead experiment the G's are even in $ while this 
is not the case for the off-diagonal G's when we deal with more than two leads. 
This conclusion, first drawn by Buttiker ^3], is seen to be quite general and not 
affected by the topology of our AB ring with two QDs. 
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5 Figure— Eight Topology 



If the two QDs are connected by a wire, the topology differs from that of a 
ring analysed so far. The neccessary modifications are quite straightforward, 
however. The AB ring is divided into two parts by the wire connecting the two 
QDs. Let $i( < i ) 2) be the flux through the lower (the upper) part, respectively. 
In a manner completely analogous to Eqs. and (|12|) . we put $i($2) onto 
W aL (W/3D respectively). With fa = 27r$ i /$ , i = 1>2, Eqs. JTTJ and |T2J) are 
thus replaced by 

Wr p = W Rp = Wp R (28) 

while 

W / LQ exp(i0i) = W / aL exp(-i(/)i) = w ha , 
W / L ^exp(-i0 2 ) = W(3 L exp(icj)2) = w L p , (29) 

with wi^p real. A further modification accounts for the presence of the wire 
which furnishes a direct link between the two QDs. We represent this link by 
a real hopping matrix element Vlr = Vrl. This element appears in the matrix 
D which now takes the form 

D PP , = S PP ,[E - Ep] + i^WppWpp, + (1 - <5ppOVr L . (30) 

p 

Except for these modifications, all formulas in Section [3] remain unchanged. 

For a discussion of the form of the matrix D in Eq. H30|) . we distinguish 
two limiting cases, where IIYrJ dominates or vice versa. It is obvious 

that for |Vrl| ^ |Tlr| we (approximately) retrieve our previous results since 
01 + 02 = 0- The distribution of the AB phase over two sets of matrix elements 
only complicates the notation. Therefore, the interesting novel limiting case is 
the one where |Vrl| ^ |Tlr|- We neglect IYr in comparison with Vrl and 
consider again the symmetric case with IY = Tr = T. Explicitly, the matrix D 
is given by 

{ v RL e-e r + (i/2)r ) ■ 

The matrix D does not depend upon the AB phase (which now appears only 
in the matrix elements W in Eq. {T3J). Moreover, the interaction Vlr causes 
standard level repulsion between the two resonances. The AB phase dependence 
of the conductance coefficients becomes complicated not because of the matrix 
D but because the electron may traverse several different paths on its way from 
the entrance channel to the exit channel. For instance, if the source (sink) is 
located in part 1 (part 2) of the AB device, there are four possible paths. One 
enters and leaves QDL, one enters and leaves QDR, one enters QDL but leaves 
QDR, and one enters QDR and leaves QDL. The relative weight of the four 
contributions depends upon the eigenvectors and eigenvalues of the matrix D. 
Again, these can be worked out straightforwardly. 
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Corrections to these limiting cases can easily be calculated in terms of a 
power-series expansion in Vlr, or in Tlr and Trl- A full diagonalization of 
the matrix D in Eq. (|3U|) is also possible, covers all the intermediary cases, and 
yields interesting results. The eigenvalues £1,2 are given by 

£1,2 = \{E h + E R - (i/2)(T L + Tr,)) 

±i ^/[e l - e r - (i/2)(r L - r R )p + 4[VL R + (*/2)t lr ][v rl + (i/2)r RL ] . 

(32) 

The eigenvalues coincide whenever the argument of the square root vanishes, 
i.e., whenever 

[E L -E r - (</2)(Tl - T R )] 2 = -4[V LR + (i/2)r LR ][VaL + (*/2)r RL . (33) 

Eq. Ij33(l extends the definition of an exceptional point to the figure eight topol- 
ogy. We note that the right-hand side of Eq. (|33J) is a periodic function of 

<j>. 

6 Summary and Discussion 

We have presented a very general approach to the transport properties of an 
AB device containing two QDs. Our main assumptions are: 

(i) The electrons do not interact. Then, we can use the Landauer-Biittiker 
approach and express every observable in terms of the conductance coefficients 
G. The latter are given as squares of the elements of the scattering matrix S. 

(ii) For the description of the the two Coulomb-blockade resonances, we use 
the single-level approximation. 

(iii) The only relevant energy dependence of S is due to the two Coulomb- 
blockade resonances, one in either QD. Then, scattering in parts 1 and 2 of the 
AB device is independent of energy and the scattering matrix S attains the 
form of Eq. ©, with [A 1 ' and independent of energy and AB phase. 

Under these assumptions, we have presented a comprehensive description of 
an AB device with the topology of a ring or of a figure eight. In particular, we 
have displayed explicitly the dependence of the S'-matrix upon energy and AB 
phase. We have shown that a novel situation arises in the case of a ring topology. 
Here the two resonances (with complex energies) are coupled via the channels 
in part 1 and part 2 of the AB device. This coupling is given by a Hermitean 
matrix multiplied by the imaginary unit i. This case differs fundamentally from 
the standard coupling of two bound states by a Hermitean interaction. The 
latter case leads to level repulsion, the former may lead to coalescence of levels. 
It seems that this phenomenon has been observed in Ref. pQ. 

We now address the approximations we have made. Perhaps most impor- 
tantly, we have neglected the Coulomb interaction between the two electrons 
populating the two QDs, and that between each of these and the electrons on 
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either QD. Inclusion of the Coulomb interaction would make it impossible to use 
the Landauer-Biittiker approach as we have done. Alternatives are discussed 
in a recent review ^j- The standard procedure employs rate equations for the 
occupation probabilities of the single-particle levels. However, this approach is 
manifestly unsuited to deal with phase correlations between scattering ampli- 
tudes. The latter are of central importance for an AB device. A more elaborate 
approach |15j uses a description in terms of an effective Hamiltonian. This 
approach assumes that the single-particle states in the QD are described by 
random-matrix theory. The effective Hamiltonian for the isolated QD is ob- 
tained as the leading term in a systematic expansion in inverse powers of g, the 
dimcnsionless conductance. The two-terminal conductance is then obtained 
from the Kubo formula and another effective Hamiltonian which includes the 
coupling to the leads. The latter is determined via a non-trivial theoretical 
derivation which in turn involves approximations. To the best of our knowl- 
edge, this approach has never been used for a multi-terminal device involving 
an AB ring. Therefore, it is not known whether the approach is able to ac- 
count for the phases which are relevant for the present system. At sufficiently 
low temperatures, the Coulomb interaction leads to Kondo-like effects in QD's. 
Remarkably, the calculation of the phase of a QD embedded in an AB ring has 
recently been worked out in the Kondo regime m spite of the difficulties 
just mentioned to deal with the Coulomb interaction outside this regime. 

In view of this situation, we can only offer a few qualitative remarks in 
support of the present approach. First, the Coulomb interaction has likewise 
been neglected in Refs. El Ej which addressed the AB phase for a single QD 
embedded in an AB ring. The results offered what seems a realistic and useful 
description of the overall phase dependence of experimental observables. Second, 
our use of the single-level approximation for each QD lends greater plausibility 
to the inclusion of the charging energy in the definition of the energies labelled 
£p. We admit, however, that the Coulomb energy between the two electrons 
(one on each QD) is not covered by this argument. Our neglect of the Coulomb 
interaction is not restricted to the neglect of the charging energy. We have 
likewise neglected the spin-dependent interaction between electrons. The latter 
is induced via the exchange term and lifts the degeneracy between singlet and 
triplet states 0ES]. This spin-dependent interaction plays a prominent role 
in Kondo-type effects. We expect that this is likewise the case in the present 
situation, especially when the two resonances overlap. Therefore, our approach 
can only be expected to work above the Kondo temperature. 

We believe that our other approximations are less severe. The single-level 
approximation should work at and near an isolated Coulomb-blockade reso- 
nance whenever resonance width and temperature are small compared to the 
charging energy. The neglect of all other energy dependence but that due to 
the resonances in the scattering matrix should be excellent barring very special 
circumstances. 
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Figure 1: Schematic representation of an AB interferometer with 5 external 
leads and two quantum dots labelled QDL and QDR embedded into its arms. 
The dotted line represents a link between the two quantum dots. 
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